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ABSTRACT 


: 1 


The effects of precession thrust pulses and energy 
dissipation upon nutation of a spin-stabi 1 iztid sj‘>acecraft are studied. 
Methods for controlling nutation during a precession maneuver are 
proposed and examined. In particular, a precession modulation control 
law is developed which uses precession thrust pulses to control 
nutation. Digital simulations show that precession control with 
separate nutation control is the fastest processing system; however, 
the precession modulation method is only fractionally slower while 
not requiring a separate nutation control system. 
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moment? of inertia aboial tlu* ] rinciral ax( s <ji 
the spacecraft, x, y, and z 

angular velocitic5; about axes f>f tlu‘ si'aeerraft 
the transverse' angular ratr 
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the angle at the start of a thrust pulse 

the time at tlie start of a thrust pulse 
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the total angle that spans those values of for 
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I. INTRODUCTION 


7\ n~stabi lized spacecraft requires a preccs>s,ion maneuver 
before reaching its final orbit. Thin maneuver in necessary to 
correct ly filign the line of action of the thruster that will be 
used for injection into the final orbit. This thruster is 
aligned with a 5 g.acecraft fixed spin axis whi cli ideally points in 
the saiiv direction as the angular momentum vector of the spinning 
spacer-raft. However, for spacecraft spinning about tin; minimum 
inertia axis, this ideal case never exists. As a result, some 
components of the total angular momentum exist along axes 
perpendicular to the spin axis. These components are characterized 
by the transverse angular rate or nutation. cince the continual 
increase of nutation will eventually cause the spacecraft to spin 
about tree maximum inertia axis, it is necessary to control this 
instai-olity by holding the nutation below some acceptable threshold 
level . 

Notation can be affected by tvjo r'ondi durino a precession 

maneuver. The first is the sequence of precession thrust pulses 
which .Ls capable of either increasing or decreasing tac amount of 
nutat.ion. The second is energy dissipation which will always 
increase the nutation of a spacecraft spinning aljout the minimum 
inorti ; axis . 

iu :ceedinq sections will discuss varir.us aspects of nutation 
control during a procession maneuver. Section II briefly describes 
the h-ar:c nutation theory. Sectioi. Ill discusses the instability 
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of nul.it i )ii during precession becriuso of energy dissipation and 
nutati'^n .'ontrol tlirough use of a separate sy.stem. Jn sr*eti(»n IV 
two procession control laws that concurrently control nutation 
are prc'sented and contrasted. Conclusions and recommendations for 
further study are given in section V. 


I 
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IT. NUTATION THEORY 

Si' i ri-st.abiliznd s^jacocraft havi- their Cartesian axi'S aJong the 

principal inertia axes with the origin at tlu? ..'enter ol m.i.ss of tlu.' 

spaceciaft. These moments of inertia are Labeled 1^, 1^ and 

are along the x, y, z axes respectively. Tlio angular velocities 

are m , , and u . Generally the angular velocity of the spin 

X y 2 

axis is considered to be constant and the equations describing the 
dynamical system are simplified accordingly. Here is considered 
the spin velocity. Thus, the angular velocities and are 
considered the transverse angular velocities and indicate the amount 
of nutation present. 

Previous nutation control studies by Grasshoff [11 and 
Taylor Ul form the basis for tLii.s paper. Tiie result ot tlisir 
studies used herein is to characterize the transverse angular rate 
or nutat-.icn, m{t), as the solution to the complex differential 
equatioTi 

w(t) + iiiio(t) = N 

where n is the nutation frequency. The applied torque- due to thrust 
pulse.s ].s related to 

N = (2) 

wnerr K is a constant and \ is the angle measured ccw from the y 

0 

axis* "\rivation <.>f oouat ion (11 is doscriLied in Appcnidix D . 



Tlie homogeneous solution to (1) is qivc^n by 


w (t) = w (t^) e 


-iitt 


( 3 ) 


Thus, U3(t) is a vector that rotates clockwise in the complex plane 
with a frequency The complete solution to (1) is given by 


^ - tj,) -11, (tj - t„)| 

,(tj) . u(t^)e ♦ ® 


(4) 


or 


■Q(t^ - t^) 


co(t^) 2 


It is noted that 


i[A - -(t^ - t^) 


:s) 




o)(to) = o>^e 


( 6 ) 


where 


u>o = w(tQ) 


(7) 


and 


“iUt^ 




( 8 ) 


^ b- 


at t 


The zero angle reference being taken as the y axis. To 
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show ui with non-zero initi^l] phase eq\iat.ion (h) 
0 


can be rewri ♦ t 'Mi as 





UIqG 


iA4. 


+ 





sin (— •) ( 



(9) 


where 


At}. = - t^) 


( 10 ) 


and 

^ = ^it . (11) 

^0 0 

The components of equation (9) are shown in Figure 1. The larger 
circle (.'-altered at the origin is the locus of the initial transverse 
veloci^.v vector. The smaller circle, or thrust lobe, tangent to the 
origin is the locus of the end of the thrust vector. The thrust 
lobe is situated as shown because of the value of a. This value is 
a constant r for all discussions to follow herein. 

The resultant vector c(t^) is a function of the variables 4)^, 
the angle at which the thrust pulse started? A4 , the duration of the 
thrust pulse; and m, the ratio of the thrust lobe radius to the 
nutation circle radius as given by the equation 


m 



( 12 ) 






Figure 1. Nufafion Circle and Prece^ion Lobe 




Figure 2. Nutation Reduction Area 


R 


In Figure 2 the interaction of these variables as tliey affect the* 
result Lo(t^) is demonstrated. The circle A with its center at l.lie 
2N ^ i(A - ^) 

tip of the vector — ^ sin(— ^)e ' describes the locus of the 

end of the vector w(tj^) as a function of 4^^, the angle at the start 
of a given thrust pulse. The interpretation of Figure 2 is that 
when the distance from the origin of the coordinate system, at the; 
center of circle B, to the circle A is less than thc‘ nutation is 

reduced. Appendix A contains a discussion of the determination of 
the angle for a given u(t^). 

This section has described nutation effects caused by a single 
thruster. The same thruster can also be used to precess a 
spin-stabilized spacecraft. For short duration thrust pulses, 
precession is given by the simple relation AH/H. 


I 
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III. NUTATION DURING AN UNMODIFIED PRECESSION MANEUVER 

In the context of this paper an unmodified procession maneuver 
is defim.-d as a maneuver in which the procession tliruster control 
signal is a sequence of pulses of duration At occurring every T 
seconds wliere T is the spin cycle pieriod. Figure 3 shows such a 
sequenrt^ of pulses. 


— »>| At j.1 T ^ 

1 1 1 1 1 1 1 .. f 

Figure 3. Precession Pulse Sequence 

As was seen in section II, thrust pulses can either increase or 
decrease the nutation of a spin stabilized spacecraft. A previous 
study by Taylor [3] shows that for zero energy dissipation nutation 
is bounded and periodic due to changes caused by precession pulses. 
The following discussion examines the bound on the nutation during 
unmodified precession with energy dissipation present. 

Equation (1) can be modified to include energy dissipation as 
follows 


( t) 


-f (- — + it) = N 

T 


(13) 


where 



JO 


T = energy dissipation time constant. (14) 

Making a notation simplif ication , (I'J) becomes 

(i) ( t ) -»■ q (i) ( t ) " N (15) 

where 


q 


1 

— + 1 ^ 

T 


(IG) 


The solution to (15) is 


-q(t - t^) -q(t - t^) 

oj(t) = ;-'(t^)e + -(1 - e ) , t < t ■■ At 

0 q 0 — — 


(17) 


, ■ ^0^ N- , -qAt -at 

a;(t) = a.'(t_)e + —(1 - e )e ' Z\t < t T 

0 q - 


( 18 ) 


with At and T equal to the duration and period of the precession 
pulses respectively. Setting t^ = 0, t = T, (18) becomes 


w[T] 


-qT N„ -qAt^ -qT 
a; 0] e + —(1 - e ) e 

q 


(19) 


or 


v.'here 


OJ [T] 


u! [0] e 


-qT 


+ F 



rq-t)e“qT 


( 20 ) 


( 21 ) 
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Now, lookJiiq at .solutions at t = 2T, t kT , 


u) [ 2T] 


u) [T] e 


-qT 


+ F 


( 22 ) 


u[kT] = (,>[ (k-DTje'*^^ 4 r (2 3) 

it is seen tliat (23) is a first order difference equation for the' 
transverse angular rate at kT as a function of the ti ansverso' 
angular rate at (k-l)T. This first order difference equation is 
solved by use of the z-transform. The z-transform of (23) is given by 

Z[w[kT]] = 2[wl (k-DTle"^'] + Z[F] (24) 

-1 -aT F 

W(z) = w[0] 4- z e W(z) 4 ^ . (25) 

1 - z 


Solving for W(z) 


W ( z ) = 


g)[0] 

-qT -1 
1 - e z 


-qT -1 

(1 - e ^ z ) (1 


z“k 


(26) 


which upon expansion becomes 


W(z) = 


D [0] 

-qT -1 
1 - e ^ z 


a - e-’-Vka 




-1 ~qT 

( 1 - z ) ( ] *- e ^ ) 


( 27 ) 


The inverse z-transform yields 
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— 1 — nnT* Fp p 

Z [W(z)] = u)[nT] = M[0]e ^ + r + 


qT -qT 

1 ~ e 1 ~ e 


(28) 


Algebraic manipulation yields 


r ^ “ qUT 1 ~ G 

:jj[nT] = u)[0]g ^ + F — 


-q(n+l)T 


1 - g 


-qT 


(29) 


or 


oj [ nT ] = [ 0 ] e 


■qnT 


+ C(1 “ e 


-q (n+1) T 


(30) 


where 


C 



(31) 


is a phasor of constant angle and magnitude. Rearranging (30) as 
follows 


:[nT] = (<o[0] - Ce e + C 


(32) 


it is evident that a: [nT] is bounded only if 


u) [0] = Ce 


-qT 


(33) 


since e is unbounded. Rewriting (32) with q replaced by (16) 


r „ ynT -iQnT 

u[nT] = Be e + C 


( 34 ) 


and 
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B = w(0] - Ce H'-.) 

clefirly shows an un.stable spiral offset from tlic origin of th< 
coordinate system by the amount of the vector addat ion ol h and C. 
This spiral is shown in Figure 4a. 

Ir. Figure 4b the nutation magnitude envelope and precej^jsion 

angle envelope are shown as a function of time. For somuj time, 

T , tlie nutation becomes so largo as to make the definition 

Critica.l 

of the direction of the spin axis meaningless for all practical 
purposi^s . 

If a precession maneuver is to be successful the nutation must 
be limited to an acceptable value. Figure 5b shows the possible 
results for a separate nutation control system output ORed with the 
precession thrust conunand signal as shown in t igure 5a. Nutation 
control systems as described in [1] and [2] would meet this 
reguirement. When a threshold value of nutation is reached thr 
system operates until the nutation is reduced suf f icieritl y . Tlie 
times when the nutation control is active are labeled , 

on the t axis of Figure 5b. A system like this will control the 
nutation build-up but has the disadvantage of complexity. It 
reguires a precession control system and an independent nutation 
control . 

The Figure 5b is only a generalized example of such a 
precess Lon-nu tat i on control sY'stem as shown in Figure 5a. He.sui ts 
of digital computer simulations fcjr a .specific spacecraft are 


shown in Appendix B. 



Figure 4a. Insfabilify of Nutation Oscillation 


precession angle envelope 


Figure 4b. Precession Angle Uncertainty and 
Instability of Nutation Oscillation 
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thruster 

command 


Figure 5a. Precession-Nutation Control Schematic 
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IV. MODULATED PRECESSION 

Section III demonstrated that a precession maneuver could be 
carried out successfully if the nutation build-up was controlled 
by a separate nutation control system. This nutation control signal 
was ORed with the precession command. In that configuration an 
independent nutation control was used to suppress the negative 
effects of both the precession thrusting sequence and the energy 
dissipation. This section demonstrates that by modulating the 
precession sequence with a gating function the precession pulses 
themselves can be used to control nutation concurrently with the 
implementation of the precession maneuver. Figure 6 shows a 
schematic for such a system. 


precession 

commands 


precession 


modulator 



thruster 

commands 


Figure 6. Precession Modulator System 


The effects of precession pulses upon nutation as functions 
of several variables are analyzed in the following. The starting 
point is equation (9) which is repeated here for convenience. 

-i(J)_ -iAi^ 2N^ * . i (X - 

u)(t-) = w_e e + - 7 ^ sin(— )e (36) 

1 U z 


I 




1 '^ 


Expandin<| into sine and cosine functions and grouping the real 

and imaginary parts separately the following result is obtained. 


2N 


(xi(t^) = u)q cos((|)q *+- A<J)) + 


0 . ,A4). ,, 

— sin (^) cos (A “ — } 


i * 


^^0 . ,A<p. . ,, 

— sin(— ^)sin(X 




- 0)^ sin((t)Q + Aif.) 


’] 


(37) 


What is sought at this point is an expression for the magnitude 
|w (tj^) j . Since 


Iw(t^) 


* 

la (t^) u) ( t^) 


(38) 


* 

where 


is the complex conjugate of cO' ' then 


(jLl(t^)l = (u) (tj^)O) (t^) ) 


(39) 


First |ui(t^) is derived by carrying out the multiplication 
indicated in (38), collecting terms, and making appropriate 
trigonometric substitutions. The result is 


I ,2 2 '*^’o . 2 A^. 

|o3(t^) I = Wq + — y- sin ( 


4N, 


0 


u)_ sin(-^) (cos(4>^ + X + ^ (40) 


The function |oj(t^) | is now known in terms of u)^ and several other 


variables . 
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Let a new function be defined as 


A(jj^ = |u)(t^)| - ]aj(tQ)| = la)(tj^)| “ 


This function describes the change in u)(t) , the nutation, due to 
precession between the times t^ and t^ . For the following discussion 
the variables A(|) and X are constants and the function 


Ao) = 


is studied. Equation (41) is now rewritten in its complete form 
which is 


% = c»o + — y- sin (— ) + “ <^0 (cos(.j)^ + X + — ) ) - 


Since Aw is to be used to control the nutation which would otherwise 
P 

increase, those values of for which 


0 > Aw = 

-- p 0 


holds true will be determined. Clearly, from (43) if the condition 
of (44) is to be met then 


0 ^ ^ sin^ (^) + w sin (^) (cos (4) + X + ^) ) . 


It should be noted at this point if 


A4^ = n27T , n = 0 , 1 , 2 , 3 , 


(46) 



19 


then no net effect on precession can be obtained. Therefore, 
those values of A4) for which 

0 < A(}) < 2v 

will be considered. Equation (45) can now be simplified to 


)sin(^) 1 cos (A + ^ 




2 - 


which becomes 


arccos 


r " o , i , . ,*♦,'1 


* 2 - 


Making a substitution based on (12) yields 


Ad) A(b . 

arccos [-m sin(— )] - X - — j ^ <|)q- 


A useful change in (50) in terms of making computations 
detailed in the following. Let 


cos ij; = -m sin(— ). 


B 0 cause of (47) the range of cos will be 

0 > cos ij/ > -1 


only 


(47) 


(48) 


(49) 


(50) 

is 


(51) 


( 52 ) 


which is shown in Figure 7. 
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Figure 7. Deferminafion of Quadrant for i/; 


Equation (52) is satisfied for 


Because of this, 


. - Stt , 

^ and — ^ ^ 


(50) can be rewritten as 


(53) 


arccos[-Tn sin (^) ] " ^ ^ 1 1 arccos[-m sin(^)] - A - 

2 II ^ u ^ III 

where arccos[-]^^ and arccos[*]^^^ refer to angles in the II and 
quadrants respectively in Figure 7. Appendix C gives graphical 
results for equation (54) as the variables m, cf)^, and A(P assvune 
various values. 

The time between the occurrence of pulses in the modulated 
precession pulse sequence determines the total amount of time 
required to complete the precession maneuver and the increase of 


2 

(54) 

III 


I 
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nutation between nutation reduction. Figure 8 provideo a graphical 
basis for a discussion of this time between modulated precession 

pulses . 


h — h- 



^ 

8Jh 







2'W 





— 1 


1 1 — 1 — Li-t 

Figure 8. Time Between Modulated Pulses 


In Figure 8 the top line represents the precession modulator gating 
function. The width of the gate is —-where <5, 

0^6 ( 55 ) 

is a function of m and as shown in Appendix C. The time for one 
nutation cycle is The sequence of vectors in the second line 

represents the leading edge of each precession pulse command coming 
into the modulator. The spacing of these vectors is T, the spin 
cycle period. For some time kT the following equation is valid. 

nT^-,(l - ^) 1 kT <_ nT,, (56 



where 
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n = 1,2,3, 
k = 1,2,3, 


The physical meaning of Figure 8 is that when a leading edge vector 
occurs in a gate window that precession pulse is used. 

Equation (56) does not permit a general analytic result for 
the values of n and k; however. Table 1 carries the results of a 
specific case. 


Gate width Maximum precession pulse 

in % of T^ separation in spin periods 


50 3 

45.5 3 

41.8 3 

36.4 3 

31.8 5 

27.3 5 

22.8 5 

18.2 5 

16.6 8 

9.1 13 

4.6 23 


Table I. Maximum Precession Pulse 

Separation for ~ 1.656 


I 
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From Table 1 it is clearly seen that as the value of 6 decreases 

the time between pulses becomes longer. 

One reason for the foregoing analysis is to determine if the 
condition can exist in which the time between gated precession 
pulses is so great that the nutation increase due to energy 
dissipation is greater than the nutation that can be removed by a 
precession pulse. A condition such as this could result in an 
unstable situation of the type discussed in section III. However, 
a long duration between precession pulses results for low 
nutation levels and, thus, stability is not degraded. 

Two methods of applying the precession modulator system are 
proposed and contrasted in the following. The first method is to 
have the precession modulator operating continuously for all values 
of 0 )^. At the start of the precession maneuver for this method of 
application the precession angle starts to increase and the 
nutation begins to decrease. As the value of decreases the 
increase in precession angle becomes less and less as a result of 
an increase in the time between allowed precession pulses. In the 
end result the increase in precession angle will become a function 
of the nutation build-up due to energy dissipation. A possible 

example is shown in Figure 9a. 

The second method of application is to use the precession 
iTjodulator only after a threshold level of is reached. Initially 
then, if the value of is below this nutation threshold the 
precession angle will increase at its maximum rate because every 
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precession pulse will be used. When the threshold value is reached 
the rate will decrease to some lesser value. A possible example 
of this method is shown in Figure 9b. 

These two methods have the following general results. The 
first method of continuous operation has the condition that the 
nutation angle is held to a very low value, but the precession 
maneuver Lakes longer. The second method of using the precession 
modulator only after a nutation threshold is reached has a faster 
rate of precession angle increase, but operates at a higher nutation 
level, a level approximately equal to the threshold value. 
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V. CONCLUSIONS 

In this research it has been demonstrated that some form of 
nutation control is required to be in operation while a precession 
maneuver is being executed. Three methods of controlling the 
nutation were examined. 

The first method consisted of a thruster command signal being 
derived from the precession control signal ORed with the output of 
a separate nutation control system. Digital simulation has shown 
this method to be the fastest of the three in attaining the final 
precession angle position. It does so, however, through increased 
system complexity and hence a possible reduction in reliability. 

The other two methods are related to each other in that they 
both use a modulated precession thrust pulse sequence to control 
the nutation. This research has demonstrated that precession 
pulses can also be used to control nutation if they are selected 
properly. This dual use of the precession pulses creates a simpler 
system because no separate nutation control system is needed. The 
first of these precession modulator methods has the modulator 
acting continuously. This method keeps the nutation level very 
low, but it also takes a much longer time to execute the precession 
maneuver. The third method which is also a precession modulator 
type, uses the modulator only after a specific nutation threshold 
is reached. This method is only fractionally slower than the 
separate nutation control system method and is still of the simpler 


I 
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precession modulator type. The nutation threshold level used in 
this third method directly affects the time required to complete 
the precession maneuver. As the nutation threshold is adjusted 
to higher acceptable levels the precession modulator acts less 
frequently and modulated precession sequence total time approaches 
that of the minimum time seen for a completely unmodulated 
precession sequence. 

In the digital simulations conducted using these methods a 
relatively small value for t, the energy dissipation time constant 
was used. This short time for T necessarily limits the value of 
the nutation threshold to a low value. 

For a precession modulator method, if the threshold level 
were set too high such that the nutation control due to the 
precession pulses could not contain the nutation build-up due to 
energy dissipation, the system would go unbounded just as if no 
attempt were being made to control the nutation build-up. Future 
study in this area of longer time constants for the energy 
dissipation can lead to an increased capability in determining 
just where the nutation threshold level can be set to avoid a 
problem such as the one mentioned above. Furthermore, the methods 
proposed above do not have any actual hardware implementation. 

The next step in that area is to generate such hardware systems as 
based upon the comprehensive theoretical analysis presented above 
which has been extensively verified through the use of digital 


simulation. 
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APPENDIX A 

Graphical Determination of 

Figure 1 of the text shows that w(t^) is a vector sum. In 
Figure 2 the circle A is the locus of the tip of vector uj(t^) as 
the precession lobe vector is held constant and the vector (.(t^) 
is allowed to assume all possible angular values (i.e. , circle A 
is the result of w(t^) as a function of . The following gives 
a graphical method for determining the angle (||q for any vector 
w(t ) from a drawing similar to Figure 2. The method is simply 
to calibrate the two circles as to angular displacement. Circle B 
is marked off in a clockwise direction (the direction of rotation 
of o)(tQ)) from the axis which is the zero angle reference. 

These angles are the values of 4 >q. The circle A is similarly 
marked off in a clockwise direction except that the zero reference 
line is shifted clockwise from a line parallel to the axis by 
an angle equal to A<J>. The angle A<(, is the angular duration of the 
precession pulse in terms of nutation circle angular measure. The 
angle for any can now be read directly from the calibration 

on circle A. If the vector u)(t^) is desired, it can then be obtained 
by drawing a radial line from the center of circle B to that point 
on circle B which is marked for the determined value of This 

graphical technique can be used to obtain the answer to equation 
(54) in the text for example. 
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An example follows in which the angles for which |w(t^)| = Wq 

are determined. The values of (J)^ for which the above condition is 
satisfied are 29° and 239° and the two u)(tQ) vector are shown 


labeled as such. 




APPENDIX B 


Digital Computer Simulation for an SMS Spacecraft Configuration 

This appendix provides the results of digital simulation of 
the nutation-precession control analysis described earlier. An SMS 
(Synchronous Meteorological Satellite) configuration with 
relatively short energy dissipation time constant is used. The 
configuration description is as follows: 

SMS Spacecraft Configuration 
Mass Properties 

Mass = 42.8 Slugs 

2 

I = 244.4 Slug-Feet 

X 

2 

I = 246.9 Slug-Feet 

y 

2 

I =97.3 Slug-Feet 

z 

o = 0.6039 

Angular Velocity 

0 ) = 90 RPM (spin axis angular velocity) 

z 

Q = 5.6916 rad/sec (nutation frequency) 

Thruster Properties 

Thrust = 5 lbs . 

Moment Arm = 3 Feet 

Duty Cycle = 1/12 seconds (4s''^ of spin period) for 
Precession 

= 0.552 seconds (1/2 nutation period) for 


Nutation 
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Results include simulation of the following nutation precession 
control laws which are displayed in the referenced figures. 

1. Figure Bl . Precession Without Nutation Control 

2. Figure B2. ORed Precession-Nutation Control 

3. Figure B3 . Full Modulated Precession (Small Initial Nutation) 

4. Figure B4 . Full Modulated Precession (Large Initial Nutation) 

5. Figure B5. Full Precession-Threshold-Modulated Precession 

The above moments of inertia are defined to be consistent with 
the description of the SMS configuration. This amounts to only a 
nomenclature change which defines as the spin axis velocity 
where as it is defined as in the text. 

In these figures the nutation level is characterized by the 
dotted line which is the plot of the actual data points generated 
in the digital simulation. Also included in each figure is the 
5.48 mr nutation angle threshold line which gives a quick reference 
for comparison of the level of nutation during each simulation. The 
5.48 mr nutation angle threshold represents the maximum amount of 
nutation that can be removed by one optimal nutation control thrust 
pulse . 

The correlation between the specific case figures of this 
appendix and comparable generalized figures in the text is good 
although exact similarity was not expected. Figure 4b is the 
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generalization of what is detailed in Figure B1 . Both show that 
as a precession maneuver is carried out without nutation control 
present the nutation increases rapidly and the direction of the 
spin axis becomes more unstable. For the case of ORed precession- 
nutation control, Figures 5b and B2 show that the precession angle 
increases steadily while the nutation level varies around the 
nutation angle threshold. Figures B3 and B4 show that the general 
Figure 9a of the text gives a moderately good description of the 
results of a continuously modulated precession. The initial value 
of the nutation is quickly brought down to a low level while the 
precession angle increase becomes very slow. Finally, it is seen 
that for the case in which the precession pulse sequence is 
modulated only after the nutation threshold is reached the slopes 
of the precession angles of Figure 9b and Figure B5 do not match 
well. It should be noted, however, that the nutation level in 
Figure B5 changes value continually which accounts for the 
numerous changes in slope of the precession angle increase in that 
figure. What Figures 9b and B5 both show is that once the nutation 
has passed the threshold level and the precession modulator starts 
acting that the slope of the precession angle becomes a function of 
the amount of nutation present. 






37 


precession •Fuel= 163.6 sec. 
Precession Angle, nutation fuel = 67. ,8 sec, 

in milliradians total fuel =271 . 4 sec. 



Figure 62, ORed Precession -Nutation Control 


Nutation Angle, in milliradians 
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Figure B4. Full Modulated Precession 
[Large Initial Nutation] 


Nutation Angle^ in milliradians 
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Precession Angle ^ 
in milliradians 


precession fuel— 226.9 sec. 
nutation fuel= 0-0 sec. 
total fuel= 228-9 sec. 
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Figure B5. Full Precession-Threshold-Modulated Precession 


I 


Nutation Angle, in milliradians 
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APPENDIX C 


Graphical Solutions to Equation (54) 


The figures in this appendix show the effects of changes in the 

variables m and A(j) upon the range of values that 4-^, the solution 

to arccosl-m sin(^)] " ^ ^ 1 <l >0 1 arccos{-m sin(^)] - A " 

2 2 U III 


A4) 

2 ' 


can take on. The value of A is held at ir radians on all figures. 

The solution for 4)^ on the figures for the above equation is the 
area(s) marked with the symbol "4 >q" • The vertical distance between 
the solution limit curves which bound 4>^ is the angle f as described 
in Figure 8 and equation (55) of the text. The conclusion of this 
appendix is that for m > 1 the duration of the thrust pulse, A4 , 


is restrictBd to less effective lengths of time. 





Figure C2. S Width m=0.6 











Figure C 




Figure Cll. 
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Development of the Complex Differential Equation 
Describing the Transverse Angular Velocity 

In this appendix a brief development of the complex differential 
equation describing the transverse angular velocity is given. 

Consider the following basic vector equation 

H = M 


where H is the total angular momentum and M is the total applied 
torque relating to a spacecraft. In terms of a Newtonian coordinate 
system and a spacecraft fixed coordinate system that rotates with 
respect to the Newtonian system the above equation is rewritten in 
terms of these coordinate systems. 


H 


N 


A + 03 x 
R — N 


(D2) 


where u) is the angular velocity between those systems. In terms of 
the total inertia I and M this becomes 


I 0 ) -f u) x(I 0 )) = M 


(D3) 


Expansion into the scalar equations gives 


lo) +(03l 03 - 03la3)=M 

XX yzz zyy X 


loj +{o3l03 - U3l03)=M 

yy zxx xzz y 


(D4) 


lO) +(o3l03 “a3l03)=M 

z 2 X y y y X X z 
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The X axis is the spin axis and o)^ = 0 is assumed. The remaining 
differential equations are rewritten as 


U) 

y 



I )0) 

X z 



03 


z 



I )03 

X y 


M 

z 

I 

z 


(D5) 


or 


where 


03 - 0) = N 

y y z y 


0) + Q 03 = N 

z Z y z 



03 



I ) 

X 


M 

N = — ^ , N 
y I z 

y 


M 

z 

I 

z 


(D6) 


(D7) 


Now by using complex varicibles notation the equations of (D6) can be 


written as 


03 (t) + iQo3 ( t) = N 


(D8) 






5^> 


where 



y 

CO (t) 

a 

= u (t) + ii^) (t) / 

V il Z 

^ z 


f 


1/2 

n = {u it ) ' , 
y z 

(D9) 


and 

a 

N = N + i(7^) N . 

V Q z 

z 



Equation (D8) is the 

same as equation (1) in the text. 





